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We introduce string-bond states, a class of states obtained by placing strings of operators on a 
lattice, which encompasses the relevant states in Quantum Information. For string-bond states, 
expectation values of local observables can be computed efficiently using Monte Carlo sampling, 
making them suitable for a variational abgorithm which extends DMRG to higher dimensional and 
irregular systems. Numerical results demonstrate the applicability of these states to the simulation 
of many-body sytems. 
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Introduction. — Explaining the properties of quantum 
many-body systems is a central topic in modern physics. 
Its difficulty is closely related to the hardness of find- 
ing a practical description of the quantum state of those 
systems. Therefore, results are usually derived by using 
either analytic approximations or numerical methods, as 
Quantum Monte Carlo (QMC) [l| or the Density Ma- 
trix Renormalization Group (DMRG) 0]. While DMRG 
works extremely well for one-dimensional systems, Monte 
Carlo proved very successful also in describing the be- 
haviour of non- frustrated systems in higher dimensions. 
Recently, an extension of DMRG to two-dimensional sys- 
tems has been developed which gives good results even 
for frustrated systems and time evolution d, 0, ■ It is 
based on Projected Entangled Pair States (PEPS) which 
describe Nature at low temperatures very well as has 
been proven by Hastings @, 0]. However, the class of 
states is too large ||, leading to an unfavourable scal- 
ing of the method in more than two dimensions or for 
periodic boundary conditions (PBC). Moreover, the al- 
gorithm relies on the underlying lattice structure so that 
irregular systems cannot be handled in a simple way. A 
subclass of PEPS which can be dealt with in an efficient 
way while keeping the power of the full family may be a 
solution to these issues. 

In this paper, we intruduce a new class of states called 
string-bond states. String-bond states form a subclass of 
PEPS which contains the relevant states in Quantum In- 
formation, as e.g. the toric code or the cluster state. Since 
expectation values can be computed easily on string-bond 
states using Monte Carlo sampling, they can be used to 
build a variational Monte Carlo algorithm for finding 
ground states. The central idea is to create the states 
by placing strings of operators on a lattice, which natu- 
rally extends the class of Matrix Product States (MPS) 
underlying DMRG to arbitrary geometries. Thus, the 
method combines the strengths of DMRG/PEPS and 
Monte Carlo: It can be applied to three-dimensional sys- 
tems, systems with periodic boundary conditions, or gen- 
eral geometries by adapting the string pattern, but also 
works for frustrated or fermionic systems which cannot 
be dealt with using Monte Carlo. At the same time, the 



computational resources scale favourably in the relevant 
parameters. We have implemented the method and suc- 
cessfully demonstrated its applicability. 

String-bond states. — Consider a classical spin system 
with configurations n = (m, . . . , tin) £ {1, . . . , d} N 
equipped with a probability distribution pin), and an effi- 
ciently computable function f(n). The expectation value 
of f(n), J2 n P( n )f( n )> can be computed using Monte 
Carlo — this is, by randomly sampling f(n) according to 
the distribution p(n) — whenever p(n) can be computed 
efficiently up to normalization. Turning towards quan- 
tum systems, for a state \tp) and an observable O 
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with p(n) — \(n\ip}\ 2 , and therefore (ip\0\ip) can be eval- 
uated using Monte Carlo whenever (n\ip) and (n\0\ip) 
can be computed efficiently. The latter reduces to (h\tjj} 
whenever O = ^ -Dfe-Ffe with diagonal and Pk per- 
mutations. In particular, this holds for local O (local 
meaning small support, as e.g. two-point correlations) 
and products of Paulis, as e.g. string order parameters. 

To build a variational Monte Carlo method, one there- 
fore has to construct states for which (n\ip) can be com- 
puted efficiently. One such class is given by Matrix 
Product States (MPS) the class of states underly- 
ing DMRG. An MPS with bond dimension D has the 
form 
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where each M*. is a D x D matrix, so that (n\ip) is 
given by the trace which can be computed efficiently. We 
generalize this to arbitrary geometries as follows. 

Definition. A state of N d-level spins is a string bond 
state if there exists a local basis \n) = |m) • • • \tin) and 
a set of strings s € S (i.e. s is an ordered subset of 
{I, . . . , TV} ) such that 
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FIG. 1: Various string patterns. In a), c) and d), the string are 
long lines - in particular, c) directly generalizes the DMRG 
ansatz, and since d) contains a), it gives a larger class of 
states. In b) and f), the strings form small loops; pattern f) 
underlies the toric code state. In e), the strings have length 
one, which suffices e.g. for the cluster state or the coherent 
version of classical Gibbs states. Patterns a), b), and both 
atop of each other have been implemented numerically. 
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FIG. 2: a) A PEPS with factorizing projectors — connected by 
the maximally entanged bonds (blue) — yields a string-bond 
state, b) To convert a general string state with many strings 
(we illustrate patterns Fig. and [lb together) into a PEPS, 
each string is routed over a separate bond. 



boundaries. Consider a PEPS with linear maps 
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for some complex D x D matrices M^ n . Here, the prod- 
uct over x £ s is over the sites x in the order in which 
they appear in the string. 

Note that the trace of a product of operators in ^ 
resembles the structure of MPS Some possible string 
arrangements are shown in Fig. [TJ 

The key point in the definition is the factorization of 
(n\ip) into efficiently computable coefficients which we 
chose to be the trace of a matrix product. However, there 
are many more natural choices, as small "blobs" with e.g. 
a PEPS parametrization or tree tensor networks (hJ. A 
special case is given by quantum states corresponding to 
thermal states of classical models [Tljj which have strings 
between neighboring sites only, cf. Fig. QJ. 

A variational ansatz based on string-bond states gen- 
eralizes DMRG beyond one-dimensional systems by com- 
bining it with Monte Carlo methods. In particular, due 
to the flexibility in the arrangement of the strings it can 
be adapted to arbitrary geometries, and the accuracy can 
be increased either by increasing D or by adding more 
strings. Clearly, the factorization of §5§ does not im- 
ply that the string-bond states themselves factorizes into 
Matrix Product States, and in fact they contain a large 
variety of relevant states. 

Properties. — Let us first clarify the relation between 
string-bond states and Projected Entangled Pair States 
(PEPS) 0]. To define a PEPS on any graph, place 
maximally entangled bonds X^c/ on eacn edge — 
associating each virtual spin with one vertex — and apply 
a linear map PM on each vertex v which maps the virtual 
spins to the o?-dimensional physical spin at v. 

For clarity, we restrict to a 2D lattice with periodic 



at site (i, j), where (<f> a \ and ((j> h \ act on two virtual spins 
each: One readily sees that together with the bonds they 
form strings (Fig.[3k), and {n\ip) is given by the product 
of the overlaps of all strings. This generalizes the states 
corresponding to classical thermal states for which P fac- 
torizes completely 11[ . On the other hand, every string- 
bond state can be written as a PEPS, even if at some 
edges many strings come to lie atop of each other. In 
that case, one places several maximally entangled bonds 
on that edge and uses one of them for each string. The 
product over the strings results in a factorizing map P 
as in ([4J, where the number of bipartite projectors {<Pi'j\ 
equals the number of strings (see Fig. [2b). Thus, the map 
pertains an efficient description, while the total bond di- 
mension scales exponentially in the number of strings. 

String-bond states are complete, i.e. every state can 
be written as a string-bond state for large enough D, 
even of form ((J) . This is easily seen by using one string 
which covers the whole system as in Fig. QJ and using 
the completeness of MPS 
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String-bond states encompass a variety of relevant 
states in Quantum Information. First, this holds for all 
MPS as e.g. the GHZ or the W state Any (general- 
ized) weighted graph state, as e.g . the cluster state [l~3j |. 
and thus all stabilizer states [14| are string-bond states 
with D = 2 [l5| . The same is true for Kitaev's toric code 
state [3], using the pattern Fig. QJ: The strings from 
loops, and the weight of each classical configuration is 
the parity of the four spins on the loop, corresponding to 
P = \Q)(ip + \(ip+\ + \l}{ip-\{ip-\ [n|. This results in a su- 
perposition of all states with an even number of |l)'s on 
any loop, which is exactly the string condensate which 
gives the toric code state. In a sense, the string-bond 
states extend the construction of Kitaev 16] and Levin 
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and Wen [ijj and may help to add new insight into Topo- 
logical Quantum Computation. Note that both the clus- 
ter and the toric code state have a block entropy which 
scales as the area, and thus string states can achieve the 
entropic area law. 

Variational ansatz. — The fact that expectation values 
of local observables and thus the energy of a local Hamil- 
tonian H can be computed efficiently on string-bond 
states allows to use them as a variational ansatz for the 
computation of ground state properties. Therefore, pick 
a string s, a site x on the string, and minimize the energy 
over the corresponding matrices {M^) n ^ =l = A (where 
A is a three-index tensor) . By iterating this protocol until 
it converges, one gets a better and better approximation 
to the ground state. 

For the optimization, we can use the linearity of string- 
bond states in A, 
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where we have explicitly denoted the dependence of 
the string-bond state \ip A ) on A. (A\X\A) denotes a 
quadratic form in A, i.e. | A) is the vectorized form of A, 
where we use boldface to avoid confusion with vectors in 
states space. Minimizing with respect to A is a gen- 
eralized eigenvalue problem and can be solved efficiently. 

However, there are a few issues which render this ap- 
proach infeasible. Firstly, X and Y have (dD 2 ) 2 degrees 
of freedom — for d = 2, D = 6, this is over 5000 — and one 
would have to do the corresponding number of Monte 
Carlo runs for a single optimization step. Moreover, 
each run would contribute to the error of X and Y", thus 
requiring a high sampling accuracy. To overcome this 
problem, we use a Monte Carlo technique called reweight- 
ing. The idea is to replace the sampling over a distribu- 
tion p{n) by the sampling over some related distribution 
p (n) ~p{n), 
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In our case, p(n) = \{n\ip A )\ 2 , Pa(n) = |(n|V>A )| 2 
(where Aq denotes the initial value of A), and f(n) — 
(n\H\ip A )/(n\ip A ), cf. Eq. (TTJ). Now define \a n ) and \b n ) 
via the linear functionals 
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Then one can readily check using (□} and © that the 
matrices X and Y in (© are 

X = J2Po(n) | b n )(a n \ , Y = Y /Po (n)\b n )(b n \ , (8) 



i.e., we can compute X and Y with a single Monte Carlo 
run. 



The second problem is the inaccuracy of X and Y 
due to the finite sampling length: In particular, errors 
in the kernel of Y will very often lead to to a wrong 
minimum. We overcome this problem by moving along 
the gradient of E(ip A ) by a small distance. With ([5]) and 
(A \Y\A ) = 1, we find that 

g^d A E^ A )\ A=Aa =Y\A ){A \X\A )-X\A ) (9) 

which only depends on absolute errors. 

At this stage, we have an applicable algorithm. An 
extra speedup of (dD 2 ) is obtained by directly sampling 
the gradient: From ©, © and (b n \A ) = 1 [Eq. Q], 



irad A E\ A=Ao 



P (n)\b n ) 



]p{m)h 



where we have defined 
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Note that for local H, h n only depends on the strings 
which intersect with H , and similarly \b n ) only depends 
on the string which contains A. As h n is independent of 
the site to be optimized, one can compute the gradients 
for all sites from the same sample, and move along all of 
them simultaneously which gives another improvement 
of the order of the lattice size. 

Numerical results. — In order to demonstrate the suit- 
ability of string-bond states for ground state calculations, 
we have implemented a simple non-opimized Matlab pro- 
gram for the string patterns Fig.QJ, (lines) and Fig.[Th+b 
(lines+loops) on a 2D lattice with PBC. Adding loops 
typically leads to a significant improvement as it gives 
full control of correlations also to the first diagonal neigh- 
bor. Additional strings which supply connections to more 
diagonal neighbors further increase the accuracy. 

We have tested our method by comparing it to the gen- 
eral PEPS algorithm 0,(3], which is the only available 
general benchmark for frustrated systems. For the frus- 
trated XX model on an 8 x 8 lattice with open bound- 
ary conditions (OBC), the general PEPS method gives 
E = —92.39 for D = 4, whereas string-bond states with 
lines+loops give E = -93.31 ± 0.02 with D = 8, while 
at the same time being about 30 times faster. Appar- 
ently, the entanglement structure of frustrated systems is 
well reproduced only for large D, suggesting that string- 
bond states are very suitable to describe such systems 
due to their favourable scaling, and since they are not 
restricted to OBC. Fig. [3] shows the magnetization for a 
frustrated XX model as a function of the transverse field 
for a 10 x 10 PBC lattice (computed with the lines+loops 
setup), for which we have no method to compare with. 

To compare the algorithm for PBC, we have there- 
fore investigated the 2D Ising model with transverse 
field and compared the results to Quantum Monte Carlo 



FIG. 3: Magnetization for a frustrated XX model (where 
each plaquette is frustrated) as a function of the transverse 
field on a 10 x 10 lattice with PBC. Note that there is no 
other method available which can deal with such systems. 

(QMC) Fig. 0] shows the magnetization and the rel- 
ative error in energy compared to QMC as a function 
of the field. Already the two basic string setups used 
reproduce both energy and magnetization very well. 

In all cases, we start with a very low number of sam- 
pling points, M = 2000, and with D = 2, and increase 
D or M or refine the gradient step adaptively. Although 
for these values M the energy is very inaccurate, the gra- 
dient is still reliable, and the method typically converges 
after about 1000 iteration steps. 

Outlook. — In this work, we have devised a class of 
states for which expectation values can be computed us- 
ing Monte Carlo sampling, and which therefore can be 
used as a variational ansatz. The central idea is that 
(n\ip) can be computed efficiently, as it can be repre- 
sented as a product of easily describable terms. In partic- 
ular, we considered string-bond states where each factor 
is a trace of a matrix product, denned on strings which 
are distributed over the system. Due to the flexibility in 
the layout of the strings, one can adapt the method to the 
geometry of the underlying system, making it an inter- 
esting approach for problems in e.g. quantum chemistry 
or irregular systems. 

The computation time scales as D 3 in the bond dimen- 
sion (D 2 for open boundaries), which improves over the 
D 5 (D 3 ) scaling of DMRG. The reason is that the tensor 
network to be contracted has dimension D rather than 
D 2 . For the same reason, Monte Carlo sampling can also 
be used to speed up the general PEPS method 0] from 
D 10 to D 6 . 

Several extensions to the ideas presented in this paper 
are being investigated. Firstly, there are many choices 
for the factors of (n\t/j) beyond strings, for instance other 
classes developed for simulating many-body systems, e.g. 
tree tensor networks [Io| which can be arranged in a way 
reflecing the geometry of the system. More generally, 



FIG. 4: Magnetization for the 10 x 10 PBC Ising model with 
transverse field: results for Quantum Monte Carlo (QMC) 
and the string setups Fig. [ij, (lines) and Fig. QJ,+b (loops). 
The right inset shows the relative error in the ground state 
energy compared to QMC as a function of the field. 

on small sets of spins ("blobs") one can allow for arbi- 
trary states or for states parametrized e.g. by a tensor 
network. Note that in all these cases, the dependence 
\ipA) on A remains linear. Finally, the string patterns 
(or general partitionings) can even depend on the classi- 
cal configuration \n), corresponding to a dependence of 
the partitionings in on the physical spin. 

Beyond the simulation of ground state properties, 
string-bond states can also be used to study time depen- 
dent phenomena or systems at finite temperature, similar 
to DMRG and the general PEPS algorithm. Finally, the 
ansatz can also be applied to fermionic systems, either 
by doing a Jordan- Wigncr transformation which does not 
affect the computability of (n\tp), or directly by using a 
fermionic ansatz for the variational Monte Carlo method. 
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Note added. — The fact that Monte Carlo can be used 
in tensor network contraction has been proposed inde- 
pendently in [2pj ] . 
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